Optical response of two-dimensional electron fluids beyond the Kohn regime: strong 
non-parabolic confinement and intense laser light 
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We investigate the linear and non-linear optical response of two-dimensional (2D) interacting 
electron fluids confined by a strong non-parabolic potential. We show that such fluids may exhibit 
higher-harmonic spectra under realistic experimental conditions. Higher harmonics arise as the elec- 
trons explore anharmonicities of the confinement potential (electron-electron interactions reduce this 
non- linear effect). This opens the possibility of controlling the optical functionality of such systems 
by engineering the confinement potential. Our results were obtained within time-dependent density- 
functional theory, employing the adiabatic local-density approximation. A classical hydrodynamical 
model is in good agreement with the quantum-mechanical results. 



Confined, two-dimensional (2D) electronic quantum 
systems have been subject to intense theoretical and ex- 
perimental investigations during the last two decadescl. 
Emerging from quantum-well structures or charge lay- 
ers in modulation doped semiconductor interfaces ,|-,they 
are nowadays routinely tailored into quantum dotsnta or 
stripsu. Possible applications range from single-electron 
transistors to coherent, tunable light sources for far- 
infrared (FIR) spectroscopy^, a method which has proven 
a powerful tool for probing slow vibrational modes in 
molecular and condensed matter systems. In the light 
of this application and in order to design future 2D THz 
devices, a detailed understanding of linear - and non- 
linear - excitation mechanisms in confined 2D electronic 
quantum systems is necessary. 

In finite 2D systems (such as, e.g., quantum dots or 
quantum strips), the linear optical response depends on 
the shape of the confinement potential vq. Recent exper- 
imental results concern parabolic or near-parabolic con- 
finement potentialso for which the so-called harmonic- 
potential theorem (HPT) states that an external dipole 
excitation can only couple to a rigid-shift mode (Kohn 
modef ) at frequency ^jK/m* independently of the ex- 
citation strcngthEl (m* is the effective electron mass and 
K the curvature of vq). In its original-formulation, the. 
HPT is a quantum-mechanical theorerroB; it was shownEI 
to hold in classical mechanics also. 

In realistic, finite 2D quantum structures, the confine- 
ment potential vo is often strongly modulated (e.g. by 
inhomogeneous charge [-distributions) exhibiting a pro- 
nounced anharmonicitylI9. Nevertheless, in maru£-.exper- 
iments, the Kohn mode dominates the response^. This 
is because for weak external fields, and for low electron 
densities, anharmonic regions of the confinement poten- 
tial are hardly explored. Experimentally it is possible to 
overcome the HPT limitations in at least two ways: either 
by increasing the density n of conduction electrons or by 
increasing the intensity of the laser light. This makes it 
possible to experimentally investigate the hydrodynamics 
of the interacting electron fluid, which is expected to re- 
veal much more information about the electron dynamics 
than the rigid-shift response in the Kohn regime. 



However, in order to adequately describe the case of 
strong fields, the non-linear response of the conduction 
electrons must b«|-epnsidered. i Wh ile the non-iijwsar re- 
sponse of atomsrTi 1 , 3 JiolcculcalJllj , nanotubestj'tZl, and 
3D quantum dotaI3'E£l has been thoroughly investigated, 
little is known about the non-linear response of finite, in- 
teracting 2D electronic systems to intense laser fields. Is 
it possible to observe higher-harmonic (HH) generation, 
either due to anharmonicities in the confinement poten- 
tial, or as a consequence of nonlinearities in the hydro- 
dynamics of the conduction electrons? Do existing THz 
sources (such as free-electron lasers) provide sufficient in- 
tensity to observe HH generation in such systems? 

Alternatively, one may consider weak fields (linear 
regime), but high conduction-electron densities. How- 
ever, most theoretical studies of the linear response ei- 
ther consider the case of low densities of conduction elec- 
trons where the_.confinement potential can be assumed 
to be parabolicEj, or the other extreme, "classical" con- 
finement by infinitely high potential barricrsEH. How the 
nature of the optical response changes in realistic sys- 
tems as the density of conduction electrons is increased 
(so that they explore more and more of the anharmonic 
parts of the confinement potential vq) is not known. How 
does the Kohn mode compete with other modes of exci- 
tation when the HPT is no longer valid? 

LasjU.but not least, to which extent can a classical 
mode£3 of the ( non-) linear response be adequate in a 
regime beyond the HPT? Such a model would have to 
account for the hydrodynamics of the electron fluid. In 
the present article we address the above questions within 
a classical hydrodynamic model, and within a quantum- 
mechapical approach, time-dependent density-functional 
theory^ (TDDFT). 

Model. Our model is described in Fig. |l|. We as- 
sume that the confinement potential vq is supplied by 
a 2D rigid positive jellium charge. If neutralised with 
charge, carriers this model corresponds to a 2D metallic 
stripEa. It can be regarded model of modulation- 
doped semiconductor heterostructures embedded in a di- 
electric medium where a layer of dopant charges corre- 
sponds to the positive background. How these, in combi- 
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FIG. 1: (a) Schematic representation of the system consid- 
ered, a 2D electron gas in the z-y-plane, confined further in 
the s-direction by a positively charged (charge density no), 
rigid jellium strip oriented along the y-axis, of width a. A 
filling fraction is defined by rj = n/no. (b) Electrostatic po- 
tential vq(x) [arb. un.] of the jellium charge (solid line) for 
a = lOag and ro ~ 0.47 [a.u.]. Reduced atomic units are used 
throughout, is the reduced Bohr radius. Around x = 0, the 
confinement potential is harmonic, vo{x) ~ (K/2) x 2 (dashed 
line) with K = 8no/a. For large values of x (\x\ > a/2), do 
grows logarithmically. 

nation with vertical gate voltages, can modify the overall 
confinement potential is discussed in Ref. [L0| 

In the following we show results for two cases, wide 
and narrow confinement in the ^-direction (a = 100 ajj 
and 10 ag, respectively); corresponding to very shallow 
confinement in the case of a = 100 a^, and very strong 
confinement for a = 10 et B . In GaAs, the widths cor- 
respond to roughly 1 /im and 100 nm, respectively. The 
filling fraction 77 is a parameter (0 < n < 1). For a given 
value of 77, the electron charge per unit length (in the 
y-direction) is taken to be the same in both cases. The 
system is subjected to an electric field E(i) = E x (t)e x 
pointing in the x-direction, e^. 

Methods. Due to translational symmetry in the y- 
direction, the problem reduces to a one-dimensional self- 
consistent one, of determining the dynamics of the elec- 
tronic density profile n(x, t) in the potential 

v([n];x,t) = xE x (t) + v KC ([n];x,t) (1) 
+ 2 Jdx'[n(x',t)-n (x')]log\x-x'\. 

Our quantum-mechanical solutionEl to this problem re- 
lies on the TDDFTEj; the exchange-correlation potential 
v^ c {[n];x, t) was treated in the adiabatic local-density ap- 
proximation (ALDA)E3. The quantum-mechanical wave- 
packet dynamics was started from the ground state of the 
unperturbed system (with E x =rQ), i- e., the solution of 
the static Kohn-Sham equations^. 

We have compared our TDDFT results with a classi- 
cal approach: the classical hydrodynamics of the electron 
fluid (r&eglecting v xc ) was solved in a co-moving Lagrange 
frameEZI, represented by a layer of infinitely thin rods of 
width da; with initial positions x, infinitesimal charges 
n(x,0)dx (per unit-length), and velocities u{x, t), evolv- 
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FIG. 2: (a) Classical (black lines) and quantum-mechanical 
(red lines) ground-state charge-density distributions for a = 
100 ajg, for different values of the parameter rj (charging frac- 
tion), from rj = 0.01 to rj = 1. (b) Same, but for a = 10 ag. 



ing according to Newton's law dtu(x, t) = —d x v([n]; x, t) 
The classical wave-packet dynamics was started from the 
stationary solution n(x) of d x v Qn]; x, 0) = 0. It was 
found by relaxing an initial Gaussian density profile with 
additional, suitably chosen Stokes damping. Since the 
static screening length of a 2D electron gas is consider- 
ably larger than the inter-particle distance, this classical 
hydrodynamical approach is expected to work well. It 
neglects exchange, correlation, and shell effects. 

Ground-state properties. Classical and quantum-me- 
chanical ground-state density profiles are shown in Fig. ||. 
For the wide system, the agreement between classical 
and quantum-mechanical profiles is satisfactory for most 
values of 77, except for large rj (where the quantum- 
mechanical profile exhibits Friedel oscillations), and for 
small values of 77 (where the quantum-mechanical pro- 
file is Gaussian while the classical one is elliptic). In the 
narrow system, the discrepancies are larger. For large 
77, electron spill-out dominates the quantum-mechanical 
profile. 

Linear response. The linear response is obtained by ap- 
plying a low-intensity white-light pulse E x (t) = Eod(t). 
A value of Eq = 0.001 [a.u.] was verified to be sufficiently 
small to remain within the regime of linear response, 
for the parameters considered here. We have calculated 
the dipolar strength function S(ui) = (2ui / Eott)Iui d(u>). 
Here d(uj) is the Fourier transform of the dipole moment. 

In Figs. ||(a) and (b) our results for S(lu) are shown. 
For low conduction-electron densities (77 = 0.1), nearly 
all dipolar strength is in the Kohn mode for both the 
wide and the narrow system, as expected. The classical 
and quantum-mechanical strength functions arc almost 
indistinguishable [insets of Figs. ||(a) and (b)]. As the 
filling fraction is increased, higher plasmon modes de- 
velop in the case of the wide system [Fig. 0(a)]. Classical 
and quantum-mechanical results agree fairly well, except 
for large values of ui where the classical plasmon disper- 
sion is found to underestimate the quantum-mechanical 
result. The results for the narrow system at 77 = 1 are 
very different [Fig_0(b)]: here we observe strong Lan- 
dau fragmentationEa of the main peak; all higher-order 
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FIG. 4: (a) Non-linear response for a = 10 a B , ljl = 
0.05 [a.u.] (0.13 THz), 77 = 0.1, and E = 3 [a.u.] (intensity 
1.38 10 6 W/cm 2 ). Shown are classical (black lines), TDDFT 
results (red lines), and quantum-mechanical independent- 
particle results (circles), (b) Density profile for the same pa- 
rameters as in (a), but Eo = 1 [a.u.] at t — 22.8 ps. (c) Time 
evolution of the width a of this density profile (the arrow in- 
dicates t = 22.8 ps). (d) Sn = n(x — x cm (t), t) — n(x) for 
a = lOOag, ljl = 0.45 [a.u.], 77 = 0.5, and Eq = 0.5 [a.u.] 
(intensity 3.8 10 4 W/cm 2 ), t = 4.5 ps. 



FIG. 3: (a) Quantum-mechanical (red lines) and classical 
(black lines) linear response for a = 100 flg and rj = 1. In- 
set: same, but for 77 = 0.1. (b) Same, but for a = 10 a B and 
?7 = 0.1, 1. (c-e) Snapshots of density profiles for a = 100 ag 
and 77 = 0.1, 0.5 and 1. For rj = 0.1, the profile is dominated 
by the Kohn mode (dashed line) . (f ) Position of the first plas- 
mon resonance as a function of rj for a = 100 03. Shown are 
classical (x) and quantum-mechanical results (□). 



plasmon modes disappear. The classical approximation 
is inadequate in this regime. 

What is the spatial profile of the modes observed in 
Figs. ||(a) and (b)? The Kohn mode (small 77) is a rigid- 
shift mode. As rj approaches unity in the wide system, 
it evolves into the first plasmon mode, a combination 
of a rigid-shift and a hydrodynamic mode (the higher 
modes are expected to be hydrodynamic modes for all 
values of rj). These two modes of oscillation correspond 
to GoldhabcpTeller and Steinwedel- Jensen modes in 
atomic nucleic^. Figs. ||(c-e) show snapshots of quantum- 
mechanical density profiles for three values of 77 in the 
wide system. By virtue of selection rules, n(x, t) — n(x) 
is antisymmetric w.r.t. reflection at x = 0. For small 77, 
the rigid-shift mode dominates and n(x, t) — n[x) oc —x 
near the origin (note that the rigid-shift profile has the 
time-dependence n(x,t) = nix — x cm (7;),0) where x cm 
is the center-of-mass of the profile). As 77 increases, hy- 
drodynamic modes emerge. They correspond, approxi- 



mately, to standing waves (reflected at ±a/2) with wave 
vectors q = -Ku/a with v = 0,1,... and frequencies 
ujp(q) = 2irn /(em*) q (e is the dielectric constant). Here 
ujp(q) is the plasmon dispersion for a spatially extended 
2D electron fluid, the confinement is modeled by assum- 
ing that the wave-length corresponding to q is given by 
the width a. Fig. ^(f) shows how the position of the first 
(y = 1) plasmon resonance evolves as a function of rj. 
For small 77, the Kohn limit is reached, as expected. As 
77 is increased, the position evolves, albeit not quite to 
the value Wp(g = 27r/a) = 27r 2 77,0 /(asm*). This is due to 
the fact that at r\ = 1, the profile of the density oscilla- 
tions is not quite sinusoidal (see also Ref. ^l[ ): sinusoidal 
modes do not diagonalise the problem, their interaction 
gives rise to a frequency shift. 

Non-linear response. The system was subjected to an 
intense monochromatic light wave of amplitude Eq and 
frequency (switched on slowly, on the time scale of 
a few cycles). The laser intensities were chosen so as to 
avoid ionisation of the system, not exceeding 10 6 W/cm 2 
(well within the range of standard free-electron lasers). 

Fig. |](a) shows classical and quantum-mechanical re- 
sults for the dipolar power spectrum |d(w)| 2 (see Ref. ||(]) 
in the narrow system. We observe excellent agreement 
between classical and quantum-mechanical results. Fur- 
ther, we observe HH at odd multiples of ojl- These HH 
are due to the electrons exploring the anharmonic po- 
tential vq (c.f. scattering of electrons off the Coulomb 
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potential in ionised atoms). The parameters (77 = 0.1 
and wl = 0.05 [a.u.]) were chosen to allow for large ex- 
cursions of n(x,t) into the anharmonic regions of vq. 

Our observations show that confined, interacting 2D 
electron fluids do exhibit HH spectra, albeit not as promi- 
nently as in single-electron systems such as atoms in 
strong laser fieldsEj. We surmise that the hydrodynamic 
modes arising from the nonlincaritics in the fluid dy- 
namics dampen the center-of-mass motion and its ac- 
celeration, reducing the intensity of HH. We have veri- 
fied that a substantial center-of-mass acceleration is ob- 
served when the non-linear electron-electron interactions 
are switched off during the laser pulse, resulting in HH 
of considerably larger strength [Fig. ||(a)]. This implies 
that independent-electron models of 2D interacting elec- 
tron fluids are likely to overestimate the strength of the 
non-linear response. Note that our quantum-mechanical 
results are not sensitive to the presence/absence of the 
exchange-correlation potential, indicating that v xc has 
little influence on HH generation. 

Following the strong external driving, the density pro- 
file n(x,t) moves with the frequency wl, but not rigidly: 
in the narrow system, the width of the profile changes pe- 
riodically (breathing mode), as shown in Figs, ^(b) and 
(c). This mode dampens the rigid-shift motion of the 
electron fluid (c.f. Ref. [u] for a similar effect in a circu- 
lar, anharmonic quantum dot). 

In the wide system (at ?; = 0.5), electron-electron in- 



teractions give rise to small-amplitude oscillations added 
to the otherwise rigidly moving density profile [Fig. ||(d)]. 
The time-dependence of the center-of-mass motion is 
found to be in good agreement with the classical model. 
In the small-amplitude oscillations, by contrast, a phase 
shift is observedo. Finally we emphasise that the selec- 
tion rules of the linear case no longer hold [Fig. f|(d)]. 

Conclusions. We have analysed the linear and non- 
linear response of confined 2D, interacting electron fluids 
to laser light. Our results may be summarised as follows: 
First, non-parabolically confined interacting 2D electron 
fluids may exhibit HH spectra under realistic experimen- 
tal conditions. HH are due to the electron fluid exploring 
anharmonicities in the confinement potential. It is found 
that electron-electron interactions dampen this effect. 
Second, with the exception of small systems at high elec- 
tron densities (where single-particle excitations interact 
with the collective modes giving rise to considerable Lan- 
dau fragmentation of the plasmons), a non-linear classi- 
cal hydrodynamical model provides a very good approx- 
imation to the linear and non-linear response obtained 
within the TDDFT; exchange-correlation effects have a 
negligible influence. It would be of interest to ascertain 
to which extent the non-linear response of geometrically 
more complex systems such as nanotubes, quantum rings, 
clusters, or Ceo can be modelled by the classical approach 
used here. 
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